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We examine the modal, spectral, and polarization entanglement properties of photon pairs gen- 
erated in a nonlinear, periodically poled, two-mode waveguide (1-D planar or 2-D circular) via 
nondegenerate spontaneous parametric down-conversion. Any of the possible degrees of freedom — 
mode number, frequency, or polarization — can be used to distinguish the down-converted photons 
while the others serve as attributes of entanglement. Distinguishing the down-converted photons 
based on their mode numbers enables us to efficiently generate spectral or polarization entangle- 
ment that is either narrowband or broadband. On the other hand, when the generated photons 
are distinguished by their frequencies in a Type-0 process, modal entanglement turns out to be 
an efficient alternative to polarization entanglement. Moreover, modal entanglement in Type-II 
down-conversion may be used to generate a doubly entangled state in frequency and polarization. 

PACS numbers: 42.65.Wi, 42.65.Lm, 42.50.Dv 
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I. INTRODUCTION 

Photon pairs generated via the nonlinear process of 
spontaneous parametric down-conversion (SPDC) can 
exhibit entanglement in multiple degrees of freedom: 
spectral, spatial, and polarization |l|. Entanglement 
arises from the multiple possibilities for satisfying energy 
and momentum conservation, as required by the para- 
metric interaction process. Spectral entanglement is ex- 
hibited equivalently in time or energy Q; spatial entan- 
glement may be manifested as entanglement in wavevec- 
tor direction transverse direction Q), orbital-angular- 
momentum (OAM) [f|, or spatial parity @, Q- A proper 
description of the quantum state of the photon pair 
(biphoton) requires a Hilbert space that includes all of 
the pertinent degrees of freedom. One of these degrees of 
freedom may be arbitrarily chosen to distinguish the two 
photons, while the other two are used to describe the 
state. The term hyper-entanglement has been used to 
describe states defined over multiple degrees of freedom 
1- 

While the polarization variable is intrinsically binary, 
spectral and spatial variables are generally continuous. 
Nevertheless, binary spectral or spatial degrees of free- 
dom can be extracted by selection of a subspacc of di- 
mension two (a qubit). Examples are frequency selection 
by use of narrow spectral filte rs p , direction selection by 
use of pinholes in the far field lfj] , and lower-order OAM- 
mode selection by suppression of higher-order modes [ill ] . 
Another approach to binarization is based on constrain- 
ing the parametric process in such a way that only two 
values of the continuous variable are permitted; such 
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configurations include narrowband spectral entanglement 
[l2| and Type-II noncollinear down-conversion in which 
only two directions (as determined by intersecting rings) 
are permitted [l3j . A third approach is to map the larger 
Hilbert space onto a binary space, as in spatial-parity en- 
tanglement 0, [HJ • 

In this paper, we consider the generation of biphotons 
by means of SPDC in a nonlinear waveguide. This con- 
figuration confines the photons to a single direction of 
propagation and discretizes the spatial degree of free- 
dom to the waveguide spatial modes. For a two-mode 
waveguide, the spatial degree of freedom is binary, rep- 
resenting a modal qubit. The waveguide also supports 
two polarizations (e.g., TE and TM in a one-dimensional 
geometry). Moreover, the spectral distributions of the 
down-converted photons may be constrained to yield a 
pair of narrow spectral lines defining a spectral qubit, or 
alternatively, may be made more flexible so as to gener- 
ate photon pairs with broad spectral entanglement; these 
conditions are achieved by making use of periodic and 
linearly chirped poling of the nonlinear medium, respec- 
tively. 

In addition to the advantage of spatial binarization, 
the two-mode waveguide has the merit of combining a 
higher rate of biphoton generation, normally obtained in 
collinear degenerate bulk SPDC, with the photon sepa- 
rability generally offered by a noncollinear configuration. 
Moreover, processing of the generated biphotons for ap- 
plications in quantum information [l5j is often facilitated 
by the use of guided- wave devices. 

In this paper, we offer a comprehensive theoretical 
study of the properties and applications of modal, spec- 
tral, and polarization entanglement of biphotons gener- 
ated via SPDC in 1-D planar and 2-D circular waveg- 
uides. In particular, we consider biphotons generated 
via different types of interactions and with the use of a 



continuous wave (CW) pump. Prior work in this area 
has been limited and has made use of a pulsed, rather 
than CW, pump 0, [13]. 

The paper is organized as follows. In Sec. II, we de- 
velop a general theory for SPDC in periodically poled 
multimode waveguides. In Sec. Ill, we determine the 
waveguide parameters required to generate modal and 
spectral or polarization entanglement in 1-D planar 
waveguides using either Type-0 or Type-II interactions. 
A similar study is carried out in Sec. IV for 2-D circu- 
lar waveguides (optical fibers). Features and applications 
of modal entanglement are considered in Sec. V. Finally, 
the conclusions are presented in Sec. VI. 



II. SPONTANEOUS PARAMETRIC 
DOWN-CONVERSION IN MULTIMODE 
WAVEGUIDES 

Consider an SPDC process in a multimode waveguide, 
where a pump wave p is down-converted into a signal 
wave s and an idler wave i . Using time-dependent per- 
turbation theory, the two-photon state \^} can be written 
as [H 

|*>~ J! drdid(r) E p (r,t) E<->(r,t) eH (r, t) |0, 0), 

(1) 

where d(r) is the second-order nonlinear-coefficient ten- 
sor; E p is the pump electric field at position r and time 

t, treated classically and assumed to be undepleted; Eq - -* 
are the negative-frequency parts of the signal and idler 
electric-field operators {q — s, i) at position r and time t; 
and |0, 0) is the vacuum state of the signal and the idler. 

For guided waves propagating along the x direction, 
the electric-field operator of the signal and idler is written 
as 



E£->(r,t) = I du q 



X exp[j(jj q t - jPm q ,a q fag) x\ 



where u q is the angular frequency; a q and m q are the po- 
larization and spatial-mode indexes, respectively; /3 m lCT 
is the propagation constant; u mqitTq is the transverse field 
profile in the y-z plane; and CT is the creation opera- 
tor for the wave q. For the classical pump electric field, a 
similar expression can be obtained by taking the complex 
conjugate of Eq. @ then replacing the creation operator 
with the wave amplitude. 

Substituting Eq. @ into Eq. (TT|), assuming that d(r) 
depends only on x, and noting that the integration over t 
yields the delta function 8 (ui p — cu s — Ui), the two-photon 
state becomes 



I*) 



dw s 



m.cr 



$m,<T (w 8 ) \cj s , m s ,a a )\u)i,mi,ai), 



where 



$m, CT (w a ) = A m ^ (w a ) / da; d(x) exp {jA/3 m>a . (u s ) x] 



A, 



{0J S 



Ay dz Yl 



(4) 

(u q ,y,z), (5) 



q— p,s,l 



m = (m s ,rrii), cr = {<7 s ,cr l ), and 



Ui = lu p - u, 

A/3 m:CT (u s ) = P mp ,a p (Wp) ~ f3m 3 ,cr 3 {u s ) ~ (3 m ,,a, (Wj) is 

the phase mismatch. The factor A mrT represents the spa- 
tial overlap of the transverse profiles of the interacting 
modes. The square-magnitude of $ mjCT (u s ) represents 
the SPDC spectrum when the signal and the idler mode 
numbers are m = (m s , rrii) and their polarization indexes 
are a = (<r s ,Oi). 

Since the generated photons are collinear, their direc- 
tions cannot be used to distinguish the two photons. If 
the mode number is used in this capacity, then the two- 
photon state described in Eq. ([3]) may be written as 



|*> ~ I duj s ^ *m, CT (W a ) |W S ,CT 



s/ m s 1^27 Ci/m 



(6) 



Alternatively, using frequency as the identifier, the two- 
photon state takes the form 



I*) 



u s ) \m s 



(7) 



We may also use polarization as the identifier and write 
a similar expression for the two-photon state. These ex- 
pressions for the state are, of course, equivalent. 

We now derive expressions for the spectral function 
&m,cr for quasi-phase matched (QPM) structures. Such 
structures are designed to phase match a specific type of 
interaction (Type or I or II) via a certain component 
of the second-order nonlinear tensor; hence d can be re- 
placed by its effective value d c g . Using electric poling, the 
nonlinear coefficient is made to alternate between ± d e g 
along one of the crystal principle axes, say x. The poling 
period can be either uniform or variable. Wc consider 
these two cases in turn. 

Uniform Poling. For uniform poling with period Ao, 
the spatial distribution of the nonlinear coefficient can 
be represented as a sum of distinct Fourier components, 

d c ft (x) = J27=i d « ex P [~0 Jo dx k k (x)] , where d K = 
(2/ttk) \d e ft\ and K K (x) = (2ttk/Aq) represent the ampli- 
tude and spatial frequency of the Kth Fourier component, 
respectively. Only that Fourier component whose phase 
is close to, or equal to, the phase mismatch Af3 m tT (oj s ) 
will contribute to <i> mi<T (u s ). Carrying out the integra- 
tion in Eq. ((4]) for a waveguide of length L, we obtain 



(w.) = A' m ^ (u s ) 



A/3 m , CT (u s )L/2tt 



(3) 



x exp j A/3 m , CT (u,) L/2 



(8) 



where 



A/3 mi0 - (uj s ) = A/3 m!<T (oj s ) 



2ir k 

^7 



(9) 



A 'm,<r i^s) = {2L/ttk) \d e s\A mta -(u s ), and sinc(6>) = 
sin(7T0) / (n9) . From a practical perspective, the pol- 
ing period is determined by satisfying the quasi-phase- 
matching condition at a certain frequency uj s , i.e., 
A/3 m , CT (w.) = 0, or A/3 m , CT (w.) = 27tk/A . 

Linearly Chirped Poling. For linearly chirped poling 
with a slowly varying spatial frequency, we can still make 
use of a Fourier series to represent the spatial distribution 
of the nonlinear coefficient with K K (x) = 2%k/Mx] = 
2ttk/Aq — £x, where C is the chirp parameter It 
is clear that it is the spatial frequency, rather than the 
spatial period, that is chirped. In this case, Eq. (Q| yields 



$, 



K) = A 'L<r exp 



-j A/3 2 (uj s 



2C 



(10) 



x [ 7 (L)- 7 (0)], 



where A'^ (uj s ) = \d eS \ A m ^ (w s ) y / 2j/n(K 2 , j(x) = 

erfi{[A0 mi<r ( w .) + (x]/V=2E}, erfi(0) = -jerf(j 0) and 
erf (•) is the error function. In the remainder of this paper, 
we use k = 1 to compute the poling periods since this 
results in the strongest nonlinear interaction and poling 
periods are not limited by fabrication techniques. 



III. MODAL ENTANGLEMENT IN 1-D 
PLANAR WAVEGUIDES 

We now consider the generation of entangled photons 
via nondegenerate SPDC in a 1-D planar waveguide. The 
structure comprises a dielectric slab of thickness h and 
refractive index n\ embedded in dielectric media with 
lower refractive index ri2, as illustrated in Fig. fT] The 
transverse profile of mode m q inside the slab is [19J 



cos 



'».„a, (Uq, V,Z) OC 



m q = odd, 



(11) 



while outside the slab it is 

Sft ( W 9> V' z ) 0Cex P(-7m a ,a 9 (Wg) z) , (12) 



where km],*, 



( fc i„ ff<! ~ #L CT J 1/2 is the z-component 



'ii... n,, — ni^q/c, 7m q ,cr q — {fi mq ,<j q 



of the wavevector k 

km a ) 1 ^ 2 j ^m q ,a q = n>2 ^>ql c i and c is the velocity of light 
in free space. Note that both n\ and are frequency- 
dependent. 

Within the slab, the mode transverse profiles are either 
even or odd functions. The propagation constants of the 
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FIG. 1: Sketch of a 1-D planar waveguide. The waveguide 
is a slab of dielectric medium of thickness h with a uniform 
refractive index n\ surrounded by media of lower refractive 
index ri2. The inner medium and the outer media are known 
as the core and cladding, respectively. The nonlinear coeffi- 
cient of the core is poled with a period A(a;); arrows indicate 
the poling direction. The poling periods at x — and x — L 
are denoted Aq and A_l, respectively. 



modes can be determined using the dispersion relation for 
the planar waveguide. For the TE wave (o-polarization), 



tan 



k'i. 



3 2 

m„,cr a 



o' 2 

II 



B 2 

In 



(13) 

whereas for the TM wave (e-polarization), the RHS of 
the above equation must be multiplied by n\/n 2 . 

Substituting Eq. (fTTj) into Eqs. dU) and ([5]), we can 
determine the functions $ m . CT (w s ), which in turn deter- 
mine the quantum state set forth in Eq. ([3]). We consider 
a two-mode waveguide for which the fundamental mode 
(m q = 0) is even and the next mode (m q — 1) is odd. If 
the down-converted photons have different spatial parity, 
the pump mode must be odd so that the spatial over- 
lap integral in Eq. ([5]), which determines A m tT (u> s ) and 
thence ^ m ,<r (w s ), does not vanish. 

We seek to generate an entangled state of the form 



|#) ~ Jduj s [$o,i )C r (Us) \w 8 ,0,a a )\ui,l,ai) 



+ $i,o,<t(w s ) \u s ,l,a s }\ui,0.,o-i)] 



(14) 



i.e., if the signal is in the even mode, the idler must 
be in the odd mode, and vice versa. This requires 
that the phase-matching condition be satisfied for each 
of these two possibilities. For a fixed geometry and 
material parameters, we find a single poling period Ao 
of the nonlinear coefficient at which these conditions, 
A/3o,i,«t (oJs) = 27r/A and A/3i :0 , CT (w a ) = 27r/A , are 
simultaneously met. We achieve this by plotting the pol- 
ing period Ao as a function of the signal frequency u> s , 
or as a function of the signal wavelength A s , for each of 
these two conditions, and search for intersections. The 
example shown in Fig. [5] reveals that this may occur at a 
single frequency or, when the two curves are tangential, 
over a broad spectral band. Once the poling period is 
selected and fixed, we determine if the spectral functions 
$o.i,ct (w s ) and $1.0,0- ( w s) overlap, i.e., have some com- 
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FIG. 2: (Color online). Values of the poling period Ao re- 
quired to satisfy quasi-phase matching as a function of the 
signal wavelength A s , for four values of the waveguide thick- 
ness h. Calculations are presented for a KTP planar waveg- 
uide with A = 0.05 pumped by a laser whose wavelength 
is Ap — 532 nm and whose modal structure coincides with 
m p = 1. The core refractive indexes are determined using 
the Sellmeier equations for KTP. The orange (black) curves 
represent condition 1 (2), namely that the signal (idler) mode 
number m 3 (m,) is and that the idler (signal) mode num- 
ber rrii (m s ) is 1. Dashed lines indicate intersection points 
or ranges, (a) Type-0 (e, e, e) with waveguide thickness h — 
1.625 /im. (b) Type-0 (e,e,e) with h = 1.49 /jm. (c) Type-II 
(o, e, o) with h = 1.617 /im. The inset shows the behavior of 
the poling period for condition 1 near the intersection point, 
(d) Type-II (e, o, o) with h = 1.56 fim. 



mon spectral band. If they do, then the two-photon state 
exhibits entanglement. 

Example. A laser source with wavelength X p = 532 
nm and in mode m p = 1, is used to pump a 1-D peri- 
odically poled KTP waveguide with fractional refractive 
index change Awl — n%/ni = 0.05. The core refractive 
indexes along the y and z axes are given by the KTP 
Sellmeier equations [20] 



n\ (y) (A) = 3.45018 + 



0.04341 



16.98825 



A 2 - 0.04597 A 2 - 39.43799' 

9 ,u 0.06206 110.80672 

< (2) (A) = 4.59423 + A2 _ - — + 

(15) 

where A is the wavelength expressed in /im. 

The curves in Fig. [5] represent uniform poling-period 



values satisfying the quasi-phase-matching conditions 
versus the signal wavelength A s , for four different values 
of the waveguide thickness h. The orange (black) curves 
represent condition 1 (2), namely that the signal (idler) 
mode number is and the idler (signal) mode number 
is 1. In panels (a) and (b), the interactions are Type-0 
(e,e,e), with |<i cff | = |d 33 | = 16.9 pm/V 0. The in- 
teractions in panels (c) and (d) are Type-II (o, e,o) and 
Type-II (e,o, o), respectively, with \d e s\ = M24I = 3.64 
pm/V. The notation (•,-,•) indicates, in consecutive or- 
der, the polarization of the down-converted photon whose 
frequency lies above the degenerate frequency, the down- 
converted photon whose frequency lies below the degen- 
erate frequency, and the pump photon. 

These plots reveal that the two curves can intersect at 
a single frequency, as in panels (a) and (c) or, alterna- 
tively, they can coincide over a broad range of frequen- 
cies, as in panels (b) and (d). Modal entanglement can 
be obtained at any intersection point since the two quasi- 
phase-matching conditions are simultaneously satisfied 
for the unique poling period determined by the point 
of intersection. In panels (a) and (c), h has been cho- 
sen such that the intersection point between the curves 
occurs near X s — 810 nm so that the down-converted 
photon with the longer wavelength lies in the telecom- 
munications window near 1550 nm. In panels (b) and 
(d) , h has been chosen such that the two curves are tan- 
gential over a broad range of wavelengths. 

We now proceed to study the spectral characteristics 
for the four structures examined in Fig. [5J The fixed pol- 
ing periods (spatial frequencies) established for panels 
(a) and (c) in Fig. [2] are used for calculating the spec- 
tra presented in panels (a) and (c) in Fig. [3J Similarly, 
the range of spatial frequencies where the curves over- 
lap in Figs. [5Jb) and (d) establish the linearly chirped 
spatial frequencies used for the calculations displayed in 
Figs. E^b) and (d). In all cases, we have chosen a struc- 
ture of length L = 25 mm. 

The output spectra for the two conditions represented 
by the quantities |$o,i.<r | 2 and I'&i^o.cr (^s) | 2 are 
shown in Fig. [31 normalized to the maximum of their 
peak values. It suffices to plot these quantities solely for 
signal frequencies above the degenerate frequency w p /2 
since 



$0, 



(16) 



This relation holds by virtue of the definitions of 
A m ,o- (w s ) and A/3 mjCT (u> s ). Modal entanglement can be 
achieved over the overlapping regions of the two spectra, 
as illustrated in Fig. [3J Linearly chirped poling leads to 
a broadband spectrum by virtue of the continuum of its 
Fourier domain, so that multiple nonlinear interactions 
can simultaneously be satisfied [22], [23[. There is, how- 
ever, a tradeoff between spectral breadth and photon flux 
density. 

The astute reader will have noticed that we have ex- 
cluded Type-I interactions from our considerations. This 
is because the down-converted photons have the same 
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FIG. 3: (Color online). Normalized output spectra 
j^m.o- (^s)| 2 = |$m i ,m i ,£r„, CTj (u s )\ 2 as a function of the sig- 
nal wavelength A s , for the four entanglement conditions 
shown in Fig. Panels (a) and (b) display modal/spectral 
entanglement regimes whereas panels (c) and (d) display 
modal/spectral entanglement and modal/polarization entan- 
glement regimes. The orange (black) curves represent con- 
dition 1 (2), namely that the signal (idler) mode number 
m s (irn) is and that the idler (signal) mode number rrii (m 3 ) 
is 1. The waveguide length L is 25 mm. (a) Type-0 (e,e, e), 
uniform poling with Ao = 7.5816 /im. (b) Type-0 (e,e,e), 
linearly chirped poling with Ao = 6.87 fim and Al ~ 6.95 
fim. (c) Type-II (o, e, o), uniform poling with Ao = 19.298 
fim. The multiple resonances that appear in condition 1 re- 
sult from the rippled behavior of the poling period near the 
intersection point, as illustrated in the inset of Fig. [21c). (d) 
Type-II (e, o, o), linearly chirped poling with Ao = 30 /im and 
Al = 70 pun. 



polarizations as those in Type-0, but Type-I suffers from 
lower efficiency. The photon flux from a Type-0 interac- 
tion is also expected to be greater than that from a Type- 
II interaction since the former exploits the strongest com- 
ponent of the second-order nonlinear tensor. 



IV. MODAL ENTANGLEMENT IN 2-D 
CIRCULAR WAVEGUIDES (SILICA FIBERS) 

Although silica is a centrosymmctric material, second- 
order nonlinearities of ~ 1 pm/V have been observed in 
poled silicate fibers, see Fig. 0] [2^ - |27| . The presence of 
a second-order nonlinearity in this centrosymmetric ma- 



FIG. 4: Sketch of the cross-section of a 2-D circular waveg- 
uide. The waveguide is made up of a cylindrical rod of dielec- 
tric medium with radius a and uniform refractive index n\ 
surrounded by a medium of lower refractive index 712 . The in- 
ner and the outer media are called the core and the cladding, 
respectively. The core is periodically poled in the x direction. 



terial appears to stem from the intrinsic third-order non- 
linear susceptibility of the glass and the built-in electric 
field arising from the displacement of charge species cre- 
ated during the poling process. Although the effective 
nonlinear coefficient |d e ff| m silica fiber is substantially 
smaller than that for many second-order nonlinear mate- 
rials (e.g., LiNbC>3 and KTP), silica fibers are nonethe- 
less of interest because they can be fabricated in very 
long lengths. Second-order interactions in poled glass 
fibers can be treated using the same techniques as those 
developed for QPM structures. 

In comparison with 1-D planar waveguides, which 
carry the index m q , each confined mode of a wave q in 
a 2-D circular waveguide is marked by an additional in- 
dex, l q . Based on the geometry of the circular waveguide, 
cylindrical coordinates are a natural choice for describing 
the transverse mode profile [l9l |: 



"i„ 



(w g ,r, <t>) 



Kz, (-yi q m q r) exp(jl q </>) 



r < a 



r > a. 
(17) 

Here r and <fi are the radial and azimuthal coordi- 
nates, respectively; l q — 0, ±1,±2, ...; m q — 1,2,3,...; 
Jz„(-) is the Bessel function of the first kind of order 
l q . Kz (•) is the modified Bessel function of the sec- 

(T) 

ond kind of order l q ; fc ; ^ 



li q m q = [Pi^-inzUq/cf] 1 / 2 ; the superscript (T) repre- 
sents the transverse component of the wavevector ki qmq ; 
a is the radius of the fiber core; and n\ and 712 are the 
refractive indexes of the core and cladding, respectively 
(they are frequency-dependent). 

The total number of modes N q is determined by the 
value of the fiber parameter V q — (uj q a/c)(nf — n^) 1 ^ 2 - 
If, for example, 2.405 < V q < 3.83, then l q = 0, ±1 and 



rn q = 1, so that N q 



3. For each value of l q , the mode 



propagation constants (3i q m q can be determined by solv- 



ing the fiber characteristic equation for k^ [fj 



X 



>la±l 



(X) 



\ (X) 



±Y- 



K 



i„±i 



(X) 



K i? (Y) 



(18) 



where X = fc( T )a and Y = {V 2 - X 2 ) 1 / 2 . The values 
of k^ that satisfy the characteristic equation are . 
Using the definitions provided above, we can then deter- 
mine Pi q m q ■ The subscript cr has been eliminated in this 
section since glass fiber is an isotropic material and the 
characteristic equation is polarization-independent. 

In 2-D circular waveguides, the two-photon quantum 
state is 



|*) ~ J du s $ tm (u s ) \uj s ,l s m s )\ui,limi), 



(19) 



where lm = (l s m s ,limi), $i m (uj s ) can be determined 
from Eq. Q and the amplitude A\ m (lu s ) in cylindrical 
coordinates is given by 

pa />27r 

Ai m (u s )= dr d(j)r TT ui q7nq (u) q , r, <p) . (20) 
Jo Jo 

Substituting Eq. JT7J) into Eq. (|20|) and performing the 
azimuthal integration leads to 5 (l p + l s + U), which can 
be viewed as a limitation on the transverse profiles of the 
interacting modes (similar to that for planar waveguides) . 

Consider nondegenerate SPDC in an optical fiber with 
2.405 < V q < 3.83. We seek to generate photon pairs 
with the following entangled state: 



(21) 



I*) ~ / du s [$oi,n K) |w„oi)Kii) 

+ $11,01 (w.) 1^,11)1^,01)], 



signifying that if the signal is in the fundamental mode 
(l s m s — 01), then the idler will be in the first mode 
{liWii = 11), and vice versa. To have nonvanishing 
A\ m (uj s ) or $i m (uj s ), the pump mode index l p m p must 
be (-11). 

Example. A pump source, with wavelength X p = 532 
nm and in mode l p m p = —11, is incident on an optical 
glass poled fiber with fractional refractive-index change 
A = 0.01. The core refractive index is computed using 
the Sellmeier equation for fused silica [l9| 



nf(A)= 1 + 



0.6962A 2 



A 2 - (0.06840) 2 
0.8975A 2 
A 2 - (9.8962) 2 ' 



0.4079A 2 
A 2 - (0.1162) 2 



(22) 



where A is the wavelength expressed in /im. Panels (a) 
and (b) in Fig. [5] represent the dependence of the uniform 
poling period on the signal wavelength A s required to sat- 
isfy quasi-phase matching for the two conditions under 
consideration, for two different values of the core radius 
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FIG. 5: (Color online). A laser source with wavelength 
A p = 532 nm and in mode l v m v = —11 is used to pump a 
poled silica fiber with A = 0.01. (a) and (b) Poling period 
versus the signal wavelength A s , with core radii a — 2.613 
and a = 2.52 /im, respectively. The orange (black) curve rep- 
resents condition 1 (2), namely that the signal (idler) mode 
number l s m 3 (hrrii) is 01 and the idler (signal) mode number 
hrrii (l s m s ) is 11. Dashed lines show the intersection point or 
region, (c) and (d) These panels display modal/spectral en- 
tanglement regimes. Normalized output spectra |$oi,n ( w s)| 2 
and |$n,oi (^s)| 2 for the structures displayed in (a) and (b), 
respectively, with L — 20 cm. The poling period of the non- 
linear coefficient is uniform in (c) with Ao = 35.813 ^m, and 
is linearly chirped in (d) with Ao = 34.9 /im and Al = 35.05 
fj>m. 



a. The curves can either intersect at a single frequency 
or they can coincide over a broad range of frequencies 
[in analogy with the results for the 1-D planar waveg- 
uide; see Figs. [U[a) and (b)]. Modal entanglement can 
be achieved where the curves intersect. The core radius 
a for panel (a) has been chosen such that the intersec- 
tion occurs at X s — 900 nm, placing the down-converted 
photon with the longer wavelength in the telecommuni- 
cations window near 1300 nm. In panel (b), a has been 
selected such that the two curves are tangential over a 
broad band of frequencies. Panels (c) and (d) in Fig. [5] 
display the normalized output spectra | $01,11 ( w s)| 2 an d 
|$n,oi (w s )| for fibers of length L = 20 cm. The poling 
periods of the fibers used in (c) and (d) are determined 
from the results obtained in (a) and (b), respectively. 
Modal entanglement can be achieved in wavelength re- 
gions where the spectra overlap. 



V. FEATURES OF MODAL ENTANGLEMENT 



In the previous sections, we have demonstrated that it 
is possible to obtain modal entanglement in 1-D planar 
and 2-D circular waveguides. In this section, we discuss 
various features and applications of modal entanglement. 

Combining the advantages of noncollinear and 
collinear-degenerate interactions in bulk crystals. 

A noncollinear configuration is often preferred in bulk 
crystals because the photon pairs generated by SPDC are 
easily separated. This process is typically not efficient, 
however, since the fraction of the photons that are en- 
tangled is seriously limited by the intersections between 
the emission cones. A degenerate-collinear configuration, 
on the other hand, creates impediments to separating the 
photon pairs but offers high efficiency by virtue of the to- 
tal overlap of the emission cones. In a two-mode waveg- 
uide with modal entanglement, in contrast, the photon 
pairs are efficiently generated in the two allowed entan- 
gled modes and are thus also readily separated by their 
mode numbers, i.e, their different transverse-field pro- 
files. This can be implemented by use of a branching 
waveguide (28[. 

Device flexibility and compatibility. Using uni- 
form or linearly chirped poling, we are able to generate ei- 
ther narrowband or broadband spectral/polarization en- 
tanglement, as illustrated in Fig. [3] By replacing the 
branching waveguide at the device output with a diffrac- 
tion grating, prism, or polarizing beam splitter (for Type- 
II), we can obtain binary modal entanglement. The 
structure can then be used as a source of binary and 
continuum entanglement that is expected to be useful 
for applications such as quantum imaging [29|-|3l|, quan- 
tum cryptography [3^1. q uantum teleportation |33|. and 
quantum information [34| . We note that this waveguide 
configuration is compatible with integrated optics, which 
can facilitate its incorporation into a practical system. 

Spectral entanglement is typically measured via 
a Hong-Ou-Mandel interferometer (HOM) (35|. An ex- 
periment is conducted by sweeping a temporal delay r 
inserted between the down-converted photons while mea- 
suring the coincidence rate of photon counts at a pair of 
detectors placed at the two output ports of the interfer- 
ometer. The coincidence rate R(t) is given by 

R(t) = ±{R -Rc{R 1 (t)}} (23) 

with 

Ro = J doj s |$ m><7 (w s )| 2 , 

Ri (r) = J du s {$ mcr <^, CT (-w.) (24) 
x exp [j (uj p - 2uj s ) t]} , 




Time delay r (ps) 



FIG. 6: (Color online). Dependence of the normalized co- 
incidence rate R on the time delay r between the two pho- 
tons calculated for a HOM interferometer, (a) 1-D planar 
waveguide with uniform poling and Type-0 interaction, (b) 
1-D planar waveguide with linearly chirped poling and Type- 
interaction, (c) 1-D planar waveguide with uniform pol- 
ing and Type-II interaction, (d) 1-D planar waveguide with 
linearly chirped poling and Type-II interaction, (e) 2-D cir- 
cular waveguide with uniform poling and Type-0 interaction, 
(f) 2-D circular waveguide with linearly chirped poling and 
Type-0 interaction. The structure and operational parame- 
ters for panels (a), (b), (c), (d), (e), and (f) are the same as 
those used for Fig. EJa), (b), (c), (d), and Fig. [5{c) and (d), 
respectively. 



where Re denotes the real part and the superscript * 
represents the complex conjugate. The integration cov- 
ers the entire spectrum, so that either ^> mo ,mi.a- (w s ) or 
3>mi,m ,cT ( w s) can be used to evaluate the integral, where 
mo and mi are the two lowest modes of a waveguide. 
Again, changing the spatial mode number is equivalent to 
changing the sign of the frequency, i.e., $ m0irahir (— u> s ) — 

The graphs presented in Fig. [5] are the normalized co- 
incidence rate expected for an HOM interferometer for 
the planar and circular waveguide structures described 
in Sees. Ill and IV, respectively. Since the signal and 
idler are nondegenerate, the familiar dip in the coinci- 
dence interferogram is modulated by an oscillatory func- 
tion. When using the mode numbers for identifying the 
photons, two different phase-matching conditions are re- 
quired to obtain spectral or polarization entanglement 
in waveguides, rather than just a single one for bulk 
crystals. This distinction arises because of the depen- 
dence of the propagation constant on the mode num- 
ber in a waveguide; in bulk crystals the magnitude of 



the wavevector is independent of the direction of prop- 
agation. Hence, A/3 ma , miiCr (uj s ) is not exactly equal 

to A0 mo>mitlJ (-w s ), so that $ m o . (uj s ) and $ m CT (-u> s ) 
do not perfectly overlap in the complex domain. This 
has several consequences: i) The interference pattern 
shifts about r = since $ mi<T (u} s ) $ m(T ( —a; s) i s a not 
purely real. This effect is illustrated in Fig. [Jjja) and 

(e) . ii) The visibility, which is defined as V (t) = 

(-Rmax — -Rmin) / (-Rmax + -Rmin), where i? max and i? m i n 

are the maximum and minimum values of the R (r), re- 
spectively, is reduced since Rq is not precisely equal to 
i? min . This is clearly observable in the Type-II curves 
portrayed in Fig. Etc) and (d). Hi) The interfer- 
ence pattern is ragged, as depicted in Fig. Hfb) and 

(f) . However, reducing the device length will result in 
a cleaner interference pattern, since the imaginary part 
of ^m.cr (w s ) CT (— u> s ) is proportional to the device 
length for both uniform and linearly chirped poling. 

Using Eqs. © and (|10[) . we define the entanglement 
length L e as the waveguide length for which the phase of 
(<^s) 'I'm CT (~ w s) is equal to 7T. For uniform poling, 
we obtain 
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FIG. 7: (Color online). Dependence of the normalized coin- 
cidence rate R on the time delay r between the two photons 
for a HOM interferometer, (a) 1-D planar waveguide with lin- 
early chirped poling and Type-II interaction, (b) 2-D circular 
waveguide with linearly chirped poling and Type-0 interac- 
tion. The structure and operational parameters for panel (a) 
are the same as those used for Figs.[3jd) and[6{d) except that 
L — 2 mm instead of 25 mm. The structure and the opera- 
tional parameters for panel (b) are the same as those used for 
Figs.[5jd) and[6£f) except that L — 20 mm instead of 20 cm. 



L K = 



2tt 



whereas for linearly chirped poling, we have 
47T 2 (A- 1 



(25) 



AZ 1 ) 



A/^ , mi , CT K)-A^ 0imiitr (-w.) 



(26) 



To explicitly illustrate that the character of the inter- 
ference pattern does depend on the device length, we 
present plots of the coincidence rate for shorter struc- 
tures in Fig. [7] The visibility of the interference pattern 
increases and its "noisiness" decreases. 

Polarization entanglement is measured by split- 
ting a pair of photons using a nonpolarizing beam splitter 
and then sending them to two Glan- Thompson analyzers 
[3o| . After passage through analyzers oriented at angles 
9\ and 62, the coincidence rate turns out to be 



Rn 



±sin 2 (29 1 ){R + Re[R 1 (T)}}) , 

(27) 

where 61 + 62 = 90°; the angles are measured 
with respect to the o-polarization wave. The visibil- 
ity of the interference pattern is given by V (r) = 
{R + Re[R 1 {t)]}/{3R - Re (r)]}. The depen- 
dence of the visibility on temporal delay is displayed in 
Figs. [8][a) and (b) for structures that support Type-II in- 
teractions. The effect of the waveguide length on the vis- 
bility is depicted in Fig.[5][c), which displays calculations 
carried out for a I-D KTP planar waveguide pumped 
using an o-polarized wave to generate Type-II (e, o, o) 
down-conversion. The nonlinear coefficient is taken to 
be uniformly poled with a period Aq = 67.7768 fim so 
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FIG. 8: (Color online), (a) and (b) Dependence of the visibil- 
ity V of the interference pattern on the time delay r between 
the down-converted photons in Type-II interactions for 1-D 
planar waveguides with uniform and linearly chirped poling, 
respectively. The structure and operational parameters for 
panels (a) and (b) are the same as those used for Fig. |3jc) 
and (d), respectively, (c) Dependence of the visibility V on 
the device length L for a 1-D planar waveguide with uniform 
poling in a Type-II interaction. The parameters are the same 
as those used for Fig. [2jd), with a uniform poling period Ao 
= 67.7768 ^m 



that the wavelength of one of the down-converted pho- 
tons lies in one of the optical telecommunications win- 
dows. The other parameters are the same as those used 
for Fig. [UJd). The visibility is seen to degrade with in- 
creasing waveguide length. Hence, there is a tradeoff 
between long waveguide length to generate a high flux of 
photon pairs and short waveguide length to generate a 
smooth interference pattern with high visibility. 



Modal entanglement can be considered as an al- 
ternative to polarization entanglement. Whereas polar- 
ization entanglement is restricted to Type-II interactions, 
modal entanglement has the merit that it can be gener- 
ated using any type of nonlinear interaction, including 
Type-0. Indeed, the down-conversion generation rate in 
Type-0 is about one order of magnitude higher than that 
in Type-II, since Type-0 utilizes the strongest compo- 
nent of the second-order nonlinear tensor. Hence, SPDC 
generated via Type-0 modal entanglement in a waveg- 
uide can serve as an efficient source of binary entangled 
photons. 

An increase in the overall signal-to-noise ratio 
in an infrared biphoton system can be achieved by 
choosing the waveguide dimensions such that one of the 
photons falls in the visible region while the other lies in 
the infrared region. This technique can be useful since 
the visible photon can be efficiently detected with a Si 
photon-counting detector while the infrared photon is de- 
tected with a less efficient InGaAs photon-counting de- 
tector. The overall signal-to-noise ratio is increased as a 
consequence (30l.[37|. 

Generation of a doubly entangled state via 
modal entanglement. Although Type-II is less effi- 
cient than Type-0, this configuration can be exploited 
to generate a doubly entangled state in frequency and 
polarization Q, 

|*) ~ + \u>i,0J s ) ® \a s ,ai) + \<Ti,a s }) . (28) 

This state can be viewed as a combination of Type-II 
(o, e, o) and Type-II (e, o, o). The photon pairs can be 
separated on the basis of their mode number with the 
help of a branching waveguide. Using Eq. ([5]), four dif- 
ferent nonlinear processes are required to generate the 
doubly entangled state: 

|*) ~ duj s [$o,i,o,e \uj s ,o) \uj i ,e) 1 

+ $l,0,o,e (w s ) \LJ i ,e) \u! s ,o) 1 (29) 
+ $o,i,e, (w s ) \u s ,e)o\uJi,o)i 

+ $1,0, e , (w s ) 1^,0)01^,6)1]. 

Although it is difficult to satisfy the phase-matching con- 
dition for all of these processes using uniform poling, this 
can be achieved by using linearly chirped poling. The 
nonlinear coefficient can also be poled in an aperiodic se- 
quence to obtain a narrowband doubly entangled state 

In Fig. IHIa), we plot the normalized output spectra 
\*&m 3 ,mi,<r a ,trt I when an o -polarized pump source in mode 
771 p = 1 at X p = 406 nm is incident on a 1-D KTP pla- 
nar waveguide with thickness h = 1.1 /im, length L = 
25 mm, and A = 0.05. The poling period of the nonlin- 
ear coefficient is linearly chirped with Ao = 6.5 /zm and 
Al = 7.9 /im. As illustrated in Fig. El the output spec- 
tra associated with Type-II (o, e, o) modal entanglement 
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FIG. 9: (Color online), (a) Output spectra of SPDC from 
a KTP planar waveguide with thickness ft = 1.1 /zm, length 
L = 25 mm, and A = 0.05, using an o -polarized source in 
mode m p — 1 at X p — 406 nm. The poling period of the 
nonlinear coefficient is linearly chirped with Ao = 6.5 /im and 
Al = 7.9 /mi. The black (1) and red (2) curves represent the 
output spectra associated with Type-II (o, e, o) modal entan- 
glement, while the orange (3) and yellow (4) curves represent 
the output spectra associated with Type-II (e, o, o) modal 
entanglement. The regions where the four spectra overlap 
indicate modal/spectral/polarization entanglement, (b) The 
regions of bulk Type-II (o, e, o) and Type-II (e, o, o) polar- 
ization entanglement are marked by small black filled circles 
and squares, respectively. The blue and orange rings repre- 
sent photons with frequencies above and below the degenerate 
frequency, respectively. 



overlap with the output spectra associated with Type-II 
(e,o, o) modal entanglement. The region of overlap be- 
tween the four spectra (except for the degenerate case) 
represents the spectral band over which a doubly entan- 
gled state can be obtained. It is worthy of note that this 
state cannot be generated in bulk crystals since four dif- 
ferent directions are required, as shown in Fig. \§lb) . In 
waveguides with modal entanglement, however, only two 
mode numbers are required. 



VI. CONCLUSION 

We have investigated nondegenerate spontaneous para- 
metric down-conversion (SPDC) in multimodc 1-D pla- 
nar and 2-D circular waveguides. Using various types 
of interactions (Type or II), the generated photons 
can be entangled in the two lowest mode numbers of 
the waveguide — we refer to this as "modal entangle- 
ment." The inherent phase mismatch in the process is 
corrected by modulating the nonlinear coefficient with 
poling periods (spatial frequencies) that are either uni- 
form or linearly chirped. Since the interaction is collinear 
within the waveguides, modal entanglement can be used 
in place of directional entanglement to distinguish be- 
tween the down-converted photons. Spectral or polariza- 
tion entanglement can therefore be generated with high 
efficiency and the entangled photons can be readily sep- 
arated. We find that there is a tradeoff between long 



waveguide length to generate a high flux of photon pairs 
and short waveguide length to generate a smooth inter- 
ference pattern with high visibility. In Type-0 interac- 
tions, frequency can be used as a mode identifier rather 
than waveguide mode number, offering the possibility of 
using modal entanglement in place of polarization entan- 
glement, which requires a Type-II interaction. This is a 
salutary feature since a Type-0 interaction exploits the 
strongest component of the second-order nonlinear ten- 
sor. The technique can be implemented in any of the 
optical telecommunications windows by controlling the 
waveguide dimensions. Finally, we have demonstrated 
that modal entanglement in a Type-II interaction can 
be used to generate a doubly entangled state in fre- 



quency and polarization. Drawing on a Hilbert space 
of higher dimensions can offer advantages in quantum- 
communication protocols and quantum computing. 
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